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FOREWORD 


,  t..  o 

The  ORO  research  group  to  which  the  authors  of  this  paper  are  assigned 

<■>  ' 

G  c  • 

just  completed  the  transportation  phase  of  a  study  of  computer-assisted  strategic 

planning  for  the  Army.  *  A  central  difficulty  in  such  planning  is  analysis  of  major 

•  o' 

transportation  networks,  so  that  one  objective  of  the  research  was  programing  a 

h '  ,  • 

digital  computer  to  carry  out  network  analysis  in  a  manner  conforming  with 

Army  requirements.  This  in  turn  called  for  the  adaptation  or  creation  of  a 

suitable  algorithm.  c 

1  The  solution  obtained  does  the  following: 

c  (a)  lor  any  existing  rail  or  highway  net,  0 
determines  the  maximum  tonnage  that  can  be 
moved  from  any  number  of  specified  "sources” 

(e.  g. ,  por+s  and  beaches)  to  any  number  of 
specified1” destinations”  (e.  g, depots  in  the 
o  Communications  Zone  of  a  military  theater  of  ,a. 

operations,  or  Army  supply  points).  ,,  o  , 

'  "  ,  Vv 

(b)  For  any  tonnage  requirements  at 
the  destinations  less  than  or  equal  to  the  max¬ 
imum,  determines  the  routings  necessary  to 

>  deliver  such  tonnage.  v'' 

(c)  Subject  to  (a)  and/or  (b)  above, 

°  determines  the  minimum- distance  routings' 

to  accomplish  the  required  deliveries. 

•-  O' 

(d)  Determines ,  for  each  routing 

in  the  solution,  The  tonnage  flow  over  sach 
individual  link.  ^  | 


♦Reported  in  ORO T- 393,  "Computer-Assisted  Strategic  Logistic 
Planning:  Transportation  Phase",  now  in  publication,  % 


o 


t,  f? 

The  computer  program  thus  obtainable  prints  routes,  with  link:  and 

junction  point  designations,  in  the  form  required  as  input  to  the  Subsequent 

■  I  'V" 

steps  of  the  transportation  planning  process  . 

A  full  discussion  of  the  application  and  some  discussion  of  the  algorithm 
is  being  made  available  in  ORO-T-393,.to  which  readers  interested  in  these  as¬ 
pects  are  referred.  It  was  decided  to  restrict  the  present  paper  entirely  to  the 

•  o 

formal  statement  of  the  mathematical  problem  underlying  the  algorithm,  with 

f  '  '  . . V  "  ;  h 

o  its  solution  and  the  required  proofs . 

The  aim  of  the  work  was  to  develop  an  algorithm  tailor-made  to  the 
Army's  requirements.  The  most  directly  applicable  ideas  in  the  literature  were 
Studied  and  adapted.  Thus  the  results  in  this  paper  are  not  entirely  original  but 
constitute  in  part  an  adaptation  of  known  methods  to  a  specific  new  application. 

J1 

Yf  <\ 

However,  the  ideal  solution*  yielding  a  minimum  "cost"  or  "distance"  flow 


pattern  for  every  quantity  of  flow  up  to  the  maximum  ,  with  the  associated  cost 
profile,  *  are  considered  to  be  novel.  ){  o 

y 

c 

k> 

I  ..  |  -  Strother  H.  Walker,  Chairman 

|  Logistics  Gaming  Group 


(n\ 


4  November  1960 

c 


♦See  definitions  below,  p  22  and  p  23  respectively. 

.tC1.  .  _  .  .  Q 
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ABSTRACT 


A  procedure  is  presented  for  solving  the  following  problem: 


,,  Given  a  finite  network  with  a  capacity 

V!  and  a  length  (or  cost  factor)  associated  with 
each  orientation  of  every  link,  find  a  family 
of  minimum-cost  flow  patterns  between  two 
specified  nodes,  one  pattern  for  every 
integral  amount  of  flow  up  to  the  maximum  ° 

consistent  With  link  capacities. 

■  \\_  ;  ' 

'  &  .  K 

C<\  » 

$)  J 

The  procedure  is  an.  iterative  process  that  adds  a  succession  of 
appropriate  elementary  chain  flows  in  such  a  way  that  each  new  pattern 

minimizes  cost  for  a  greater  amount  of  flow.  The  presentation  is 

« . *  Q  ■  °  ■, 

intended  to  be  self-contained  and  includes  a  proof  of  the  validity  of  the 

procedure. 
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A  PROCEDURE  FOR  DETERMINING  A  FAMILY  OF 
MINIMUM-COST  NETWORK  FLOW  PATTERNS 


O 


INTRODUCTION 

■■  c 

This  paper  presents  a  general  procedure  for  solving  the  following, 
abstract  network  flow  problem:  Given  a  finite  network,  whose  links  have-as¬ 
signed  capacities  and  lengths  (Or  cost  factors)  ,  find  a  family  of  minimum  cost 
flow  patterns  between  two  specified  nodes,  one  pattern  for  every  integral 
amount  of  flow  up  to  the  maximum  consistent  with  link  capacities. 

The  first  section  defines  the  basic  network  terms  and  symbols  employed 
throughout  the  paper,  A  number  of  lemmas  that  are  consequences  of  these 
definitions  are  also  established  here.  The  second  section  presents  an  iterative 
procedure  that  is  shown  to  solve  the  stated  network  problem,  under  the  as¬ 
sumption  that  an  algorithm  with  certain  characteristics  can  be  devised.  The 
final  section  of  the  paper  produces  this  algorithm  and  proves  that  it  has  the 
required  characteristics.  A  glossary  of  the  principal  symbols  appearing  in 
the  development  is  included  at  the  end. 

Some  of  the  basic  ideas  introduced  in  the  paper  are  illustrated 

geometrically.  However,  the  mathematical  development  does  not  rest  on 

>)  ~  ... 

geometrical  or  topological  considerations.  It  should  also  be  noted  that  net--  / 
work  theory,  including  the  specific  topic  of  this  paper,  is  applicable  to 
problems  in  which  the  relevant  network  is  not  a  physical  structure  connecting 

.  ii  il 

points  in  space.  Whereas  network  links  may  be  interpreted  as  transport  routes, 
they  may  also  be  interpreted  as  the  set  of  permissible  information  channels  in 
an  organization  or  as  the  set  of  possible  transitions  between  the  various  states 

f  i 

of  a  system. 


PRELIMINARY  DEFINITIONS  AND  LEMMAS 

s 

BASIC  NETWORK  TERMS 

A  "network”  consists  of  a  finite  set  of  distinct  elements  N. ,  N0,  .  .  . ,  N, 

<  \\  . 1  £ k 

called  "nodes"  ,  together  with  a  subset  of  all  unordered  pairs  (N. ,  N.)  that  can  be 
formed  from  distinct  nodes  N.  and  N...  The  elements  of  this  subset  are  termed 
the  "links"  of  the  network.  These  are  necessarily  finite  in  number,  an  upper 
limit  on  the  number  of  links  being  l  k(k-l)]/2,  where  k  is  the  number  of  nodes. 
(This  is  the  number  of  links  if  every  pair  of  distinct  nodes  determines  a  link.  ) 

If  an  ordering  is  assigned  to  the  nodes  of  a  link  (N.,  N.)  designating  N.  as 

•>  ,  ..  1  J  J 

the  first  or  "ihitiaLnode"  and  N.  as  the  second  or  "terminal  node",  .  ,  the  re- 

0  C-  1 

suiting  object  is  termed  a  "directed  link  from  Ny  to  N.  and  is  denoted  by  N.;,  N. . 

.J'  ^  3  ^ 

Thus  two  directed  links;,  N.,  N.  and  N. ,  N. ,  can  be  associated  with  each  link. 

It  is  convenient  to  introduce  a  representative  set  of  directed  links  called 
an  "enumerating  set.”  If  a  network  has  n  undirected  links,  an  enumerating  set 
for  this  network  is  any  set  of  n  directed  links  such  that  each  undirected  link  is 
represented  by  precisely  one  of  the  two  directed  links  associated  with  it.  It  is 
evident  that  there  are  211  distinct  enumerating  sets  that  can  be  associated  with 
a  network.  The  function  of  such  a  set  is  to  establish  for  each  link  a  convention 
as  to  which  direction  will  be  considered  positive  in  describing  flows  through 

■  :  --  "  ;i. 

the  network.  -  -  c 

b  Any  network  can  be  represented  by  a  simple  geometrical  structure  in 

three-dimensional  space.  Let  any  k  distinct  points  represent  the  k  nodes  .  For  ^ 


each  link  (N. ,  N.)  of  the  network  construct  a  simple  curve,  having  the  corresponding 
points  as  end  points,  that  contains  none  of  the  k- 2  points  corresponding  to  other 
nodes  and  that  intersects  none  of  the  curves  representing  other  links.  Such  a 
structure  is  termed  a  "linear  graph/”  Figure  la  illustrates  a  simple  network 
having  6  nodes  and  8  links;..  C;  ([ 

If  A  and  B  are  distinct  nodes  of  a  network,  an  "A/B  path"  is  a  finite  set 

of  directed  links  that  can  be  assigned  an  ordering  M.,  N.,  i  =  1,  2,  .  ,  .  ,  m  (m 

<> 

being  the  number  of  links)  so  that  =  A,  Nm  =  B  and  N.  =  M.+^  for 


i=  1,  2,  ,  .  .  ,  m-1.  A  and  B  are  termed  the  "initial"  and  "terminal"  nodes 
of  the  path  respectively.  Figure  lb  represents  a  path  from  Ng  to  Ng  since  the 


directed  links  can  be  arranged  1:  the  following  sequence;  Ng,  Ng,  N^,  N2,  N2,  N^_, 


N 


^,.N^,  N^.Ng,  N2,  Ng.  A  network  is  said  to  be  "connected"  if  for  every  pair  Of 
distinct  nodes  A  and  B  there  exists  at  least  one  A/B  path.  Note  that  if  a  network 


is  connected  in  this  sense  a  linear  graph  Constructed  to  represent  the  network  is 


a  connected  point  set  in  the  conventional  sense  of  connectedness ,  In  the  remainder 
of  this  paper  all  networks  under  consideration  are  assumed  to  be  connected. 


A  "loop"  in  a  network  is  a  finite  set  of  directed  links  that  can  be  assigned 


an  ordering  M.,  N.,  i  =  1,  2,  .  ,  .  ,  in,  such  that  M1  =  Nm  and  N.  -  M.+1  for 
1=1,  2,  ...  ,  ,  m-1.*  Figure  lc  is  an  example  of  a  loop.  .. 

Note  that  the  same  directed  link  may  appear  more  than  once  in  an  A/B 
path  or  a  loop  and  that  both  M,  N  and  N,  M  may  appear.  If  the  same  directed 

link  appears  k  times,  it  is  considered  as  being  enumerated  k  times  in  the  set 

•  -  "  : 

of  links  that  defines  the  path  or  loop. 

C  -  ■  ■'  o 

♦Note  that  a  loop  could  be  considered  as  an  A/A  path.  In  this  paper  it 
is  never  categorized  as  such,  however. 


5 


m 


d.  Ng/Ng  chain  «  e.  A 


Fig.  1-— Illustration  of  Basic  Concepts 


<> 


Since  the  links  that  form  a  path  or  loop  may  be  arranged  in  a  sequence  of 
the  form  M^Mg,  M^M^,  .  .  .  ,  it  is  frequently  convenient  to  represent 

a  path  or  loop  by  the  more  concise  notation  M^,  Mg,  ,  .  ,  M^,  where  the  notation 

represents  an  M^/M^  path  if  4  and  a  loop  if  -  M^.  It  should  be  noted 
that  this  representation  is  not  in  general  unique,  however.  For  example 
Mv  Mg,  Mg,  M4>  M2,  Mg  and  M^  M4>  Mg,  Mg,  M^,  Mg,  Mg  both  represent  the 
same  M./M_  path,  since  they  are  composed  of  the  same  6  directed  links. 


15! 

An  "A/B  chain'’  is  defined  to  be  an  A/B  path  such  that  the  initial  nodes 

(and  consequently  also  the  terminal  nodes)  of  its  directed  links  are  all  distinct. 

r  <>  ^ 

Similfl/.y  a  '’cycle"  is  a  loop  such  that  the  initial  (or  terminal)  nodes  of  its 

links  are  all  distinct.  If  interpreted  geometrically,  chains  and  cycles  correspond 
to  simple  or  nonintersecting  curves  in  a  linear  graph.  Figures  Id  and  le 

illustrate  an  N„/N_  chain  and  a  cycle  respectively.  As  a  result  of  these  defi- 

v  a  <n, 

nitions,  a  chain  has  a  unique  representation  of  the  form  M^  Mg,  .  .  .  ,:  M^, 
and  the  representation  of  a  cycle  is  Unique  except  for  cyclic  permutation  of  its 
links.  (For  example,  M1,Mg,  Mg,M^,M1  and  Mg,  M^.  M^  Mg,  Mg  represent 
the  same  cycle. ) 


CAPACITATED  AND  WEIGHTED  NETWORKS 


If  a  network  has  k  undirected  links,  and  a  nonnegative  integer  termed  a 
"capacity"  is  associated  with  each  of  the  2k  directed  links,  the  network  is  termed 
a  "capacitated  network."  The  capacity  of  a  directed  link  M,N  is  denoted  by 
c(M,  N).  The  capacities  need  not  be  symmetric,  i.  e. ,  it  is  permissible  that 
c(M,  N)  4  c(N,  M).  The  link  capacities  serve  as  upper  limits  on  the  rate  of  flow 


through  the  corresponding  links  in  the  subsequent  discussion.  If  the  network  is 

a  mathematical  model  of  a  transportation  network, . .  capacities  might  have 

the  dimensionality  ’’tons  per  unit  time."  '• 

If  a  nonnegative  real  number  called  a  "length"  is  associated  with  each  of 
the  2k  directed  links  of  a  network,  the  network  is  considered  "weighted. "  The 
length  of  a  directed  link  M,  N  is  denoted  by  1(M,  N).  Here  again  the  lengths  need 
not  be  symmetric.  These  lengths  serve  as  Weighting  factors  applied  to  the 
amount  of  flow  in  each  link  in  the  subsequent  discussion.  Depending  on  the  ap¬ 
plication,  the  length  of  a  directed  link  may  be  an  abstraction  of  physical  length, 
cost  per  unit  of  flow  from  the  initial  to  the  terminal  node  of  the  link,  or  some 
other  attribute . 

//  ‘  -N  %  0* 

<■  * 

If  a  represents  a  path  or  loop  in  a  weighted  network,  the  "length  of  a", 

c  , , 

denoted  by  1(a),  is  defined  as* 

m  =  ^i(mTn) 

a 

In  particular  if  a  is  the  W  N  path  consisting  of  the  single  directed  link  M,  N, 
1(a)  =  I(M,  N);  thus  there  is  no  numerical  ambiguity  in  using  the  symbol  "1" 
for  both  the  concepts  of  link  length  and  path  or  loop  length. 

,c  “  * 

O  ; 

*In  this  equation  and  in  other  places  where  there  is  no  ambiguity,  the 
notation  S  is  used  in  place  of  the  more  cumbersome  notation  S  . 
a  M,  Nea 


o 


LINK  FLOWS  AND  FLOW  PATTERNS 


Consider  now  the  concepts  of  link  flows  and  flow  patterns  in  a  network, 

A  flow  in  a  link  is  a  nonnegative  integer  p  associated  with  a  particular  orienta¬ 
tion  of  the  link.  If  the  flow  is  designated  as  being  from  M  to  N  in  link  (M,  N),  a 
flow  f  is  designed  for  each  of  the  directed  links  M,  N  and  N,  M  as  follows: 

LS  f(M7N)  =  p;  f(W)  =  -p 

Thus  the  flow  value  associated  with  a  directed  link  is  negative  if  the  direction  of 
the  link  is  counter  to  the  assigned  direction  of  the  flow;  otherwise  it  is  positive. 

_ _  O 

Note  that  since  f(M,  N)  =  -f(N,  M)  it  is  sufficient  to  prescribe  the  values  of  f  for 
an  enumerating  set  in  order  to  establish  flows  for  all  directed  links  of  the  net- 
work.  If  f(M,  N)  =  f(N,  M)  =  0,  it  is  still  Convenient  in  some  contexts  to 
arbitrarily  consider  the  flow  as  being  "from”  one  node  "to”  the  other. 

If  a  flow  is  assigned  to  a  link  (M,  N)  of  a  capacitated  network,  the  flow 


is  said  to  be  "feasible"  if 


-c(N,  M)  <  f(M,  N)  <  c(M,  N) 


or  equivalent  if 


-c(M,  N)  <  f(N,  M)  <  c(N,  M) 


If  f(M,  N)  -  c(M,  N),  the  directed  link  M,N  is  said  to  be  "saturated"  by  the  flow. 

A  flow  saturates  both  directed  links  M,  N  and  N,  M  only  in  the  trivial  case  where 

c(M,  N)  =  c(n7M)  =  f(M,  N)  =  0* 

♦For  if  f(MTN)  =  c(M,  N)  and  f (N,  M)  =  c(N,M),  it  follows  that  c(M,  N)  +  c(N,  M)  =  0, 
since  f(M,  N)  +  f(N,  M)  =  0.  But  since  both  capacities  are  nonnegative,  both  must  be  zero. 


/> 


o 


It  is  important  to  stress  the  fact  that  one  cannot  independently  assign 
flows  to  each  directed  link  of  a  pair,  although  it  is  evident  from  the  preceding 
definitions,  because  f(M,N)  and  f(N,M)  are  two  equivalent  representations  of 
the  same  flow,  which  is  best  thought  of  as  a  variable  assigned  to  the  undirected 
link  (M,N). 

Consider  a  network  in  which  a  flow  is  assigned  to  every  link.  Then  if 
a  node  M  is  fixed  and  any  enumerating  set  S  is  selected,  the  expression 


M  fixed 


represents  the  total  flow  away  from  node  M  minus  the  total  flow  into  node  M  for 

\  ^  , 

those  links  whose  representative  in  S  has  M  as  initial  node.  Similarly 


^f(P,M), 


M  fixed 


represents  the  total  flow  into  node  M  minus  the  total  flow  away  from  node  M 
for  those  links  whose  representative  in  S  has  M  as  terminal  node,.  Thus  the 
function 


Y(M)  =2^  f(M,N)  -  ^  f(P,M),  M  fixed 


represents  the  total  flow  away  from  node  M  minus  the  total  flow  into  node  M. 
Y(M)  is  termed  the  "net  output  at  node  M."  The  function  Y  is  independent  of 
the  enumerating  set  S  since  f(M,N)  =  -f(N,M).  If  the  representation  in  S  of  a 


10 


link  is  changed  from  M,  N  to  N,  M,  the  effect  is  merely  to  transfer  a  term  from 
one  of  the  above  summations  to  the  other,  with  a  corresponding  change  of  sign. 

If  a  flow  is  assigned  to  every,  link  of  a  network  and  Y(P)  =0  at  every  node 
P  except  possibly  A  and.'7?,,  then  it  is  readily  shown  that  Y(A)  =  - Y(B) ,  because 


^f(M,  N)  -  Yf(N,  M)  =  0 


since  each  M,  NcS  appears  once  in  each  summation.  Also  SY(P)  -  Y(A) ,  which 

Pt=B 

can  be  written  as 


^f(M,  N)  -  ^f(N,  M)  =  Y(A) 


S 

M1=B 

Subtracting  Eq.  2  from  Eq. 


S 

Mt-B 


Jf(M,N)  -  Jf(N,M) =  -Y(A) 


S 

M=B 


•  \  S 

0  v  M=B 


but  the  lefthand  side  is  Y(B)  by  definition.  If  Y(A)  >  0,  the  assignment  of  flows 
is  termed  an  "A/B  flow  pattern. "  Y(A)  is  called  the  '’value”  of  the  flow  pattern, 
A  is  termed  the  "source",  and  B  is  termed  the  "sink. "  If  F  represents  an  A/B 
flow  pattern,  the  symbol  V(F)  is  also  used  to  denote  its  value.  V(F)  is  neces¬ 
sarily  an  integer  since  all  link  flows  are  integers  by  definition.  If  Y(P)  =  0  for 
every  node  P,  the  assignment  of  flows  is  termed  a  "zero  flow  pattern. "  (In 

o'1  ■ 1  ' 

some  contexts  it  is  convenient  to  consider  a  zero  flow  pattern  as  an  A/B  flow 

Jj  ,  O  \\  ^  ^ 

pattern  whose  value  is  zero. )  To  summarize,  an  assignment  of  flows  to  the 


V. 


links  of  a  network  is  a  flow  pattern  in  this  paper  if,  and  only  if,  Y(P)  0  for  at 


most  two  nodes. 

If  P. ,  F„,  .  .  .  ,  F  are  A/B  flow  patterns  (some  or  all  of  which  may 
have  the  value  zero)  and  if  ^(M,  N),  f2(M,  N),  .  ......  ,  fn(M,  N)  represent  the 

corresponding  flows  in  an  arbitrary  directed  link  M,  N,  then  the  patterns  are 


said  to  be  ’’conformal”  if  for  every  link  M,  N  either  fv(KT,  N)  >  0  for 
i  -  i,  2,  .  .  ...  ,  nor  f.(M,N)  <0  for  i  =  1,  2,  .  ....  .  ,  n.  Expressed  differently 


they  are  conformal  if  there  is  no  link  such  that  there  is  nonzero  flow  in  one 


direction  relative  to  some  F.  and  nonzero  flow  in  the  other  direction  relative  to 


some  other  F 


r 


An  algebra  of  flow  patterns  can  be  defined.  If  F  and  G  are  A/B  flow 

'  'l, 

patterns  (where  the  possibility  of  zero  flow  patterns  is  included)  and  f  (M,  N)  and 
g(M,  N)  denote  the  corresponding  flows  in  an  arbitrary  directed  link  M,  N,  the 


sum  of  F  and  G,  denoted  by  F©  G,  is  defined  as  follows;  for  every  directed 

___  0  _ _  J 

link  M,  N  let  h(M,  N)  =  f(M,  N)  +  g(M,  N).  The  values  of  h  then  determine  the 

)} 

link  flows  of  F©  G.  Clearly  F©  G  is  also  an  A/B  flow  pattern,  whose  value 
is  V(F)  +  V(G). 

Similarly  F0G  is  defined  by  the  relation  h(M,  N)  —  f(M,  N)  -  g(M~N)  .  In 
this  case  if  V(F)  >  V(G),  F©G  is  an  A/B  flow  pattern  of  value  V(F)  -  V(G).  In 
the  contrary  case  F0  G  is  a  B/A  flow  pattern  of  value  V(G)  -  V(F). 

Elementary  Flow  Patterns  °  c. 

//  ’  ''  '  ft . 

Certain  elementary  flow  patterns  corresponding  to  paths  and  loops  play  a 
central  role  in  the  subsequent  development.  If  P  is  an  A/B  path  or  loop,  if  a 


'// 

-r/ 

;  V  .  .  0- 

♦  — - - 

positive  flow  of  m  units  is  assigned  to  each  link  M,N  of  P  (  and  necessarily  also 

a  flow  of  -m  units  is  assigned  to  N,  M),  and  if  zero  flow  is  assigned  to  all  other 

links  of  the  network,  then  this  assignment  is  called  an  ”A/B  path  (or  loop)  flow” 

of  m  units.  If  m  =  1,  it  will  also  be  designated  a  ’’unit  path  (or  loop)  flow.”  The 

.  ,  '  /  r  ’Z 

notation  m(P)  Or  m(M^,  .  .  .  ,  M^)  is  used  to  denote  such  an  assignment  of 

flows.  If  P  is  an  A/B  path,  it  is  readily  seen  that  m(P)  is  an  A/B  flow  pattern 
of  value  m.  If  P  is  a  loop,  m(P)  is  a  zero  flow  pattern.  Chain  and  cycle  flows 
are,  of  course,  special  instances  of  path  and  loop  flows  respectively.  The  fol¬ 
lowing  lemma  demonstrates  that  a  unit  A/B  path  flow  can  always  be  decomposed, 

)\ 

‘in  a  certain  sense,  into  simpler  components . 

c>  .O' 

Lemma  1:  If  P  is  an  A/B  path,  then  1(P)  =  1(C)©  Z* 

where  C  is  an  appropriate  A/B  chain,  Z  is 
a  zero  flow  pattern,  and  1(C)  and  Z  are 

conformal.  „  o 

_ -  _ _ .  .  ^ 

Proof:  P  can  be  represented  as  N  N0,  .  .  .  ,  N.  .  where  N.  =  A  and 

1  &  K  1  vV/ 

~  B.  If  P  is  a  chain  the  result  follows  immediately  by  taking  C  =  P  and 

o 

letting  Z  be  the  flow  pattern  that  consists  of  zero  flow  in  every  link  of  the  net¬ 
work.  If  P  is  not  a  chain,  there  is  a  first  index  i  such  that  N.  ==  N.  for  some 

<  ■  c  '■ .  1  l 

j  <  i.  It  is  readily  seen  that  1(P)  =  1(P1)©1(S,)  where  P.  is  the  A/B  path 

I  1  1  )■<  ^ 

N^7  •  .  ,  ,  Nj7  N.+1,  .  .  ,  ,  Nk  and  is  the  cycle  K,  N.+  y  .  .  ,  N.. 

♦If  a  given  link  M,  N  appears  more  than  once,  say  R  times,  then  m*  R 
units  of  flow  are  to  be  asssigned  to  M,  N.  \  ^ 
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If  P_1  is  a  chain,  the  proof  is  completed  hy  setting  C  ~  and  Z  =  If  P. 

is  not  a  chain,  the  above  process  can  be  repeated,  expressing  i(P1)  as  l(p2)©l(S2) 
so  that  i(P)  =  1(P_)©  1(S1  )03-(So).  So  long  as  each  successive  path  P  obtained 
by  this  process  is  not  a  chain,  it  is  always  possible  to  "extract"  another  path 

pmfi eo  that  i<v  ■  ‘iwenyw  KP)  -  .  •  @i«sm+1) 

% 

But  clearly  the  process  must  terminate  since  P^+ ^  has  fewer  links  than  P  ( .  It 

follows  that  at  some  stage  M,  PM  is  a  chain.  Then  if  C  =  PM  and 

Z  ?=  1($1)0  .  .  .01(SM),  one  has  the  desired  expression  for  1(P).  Z  and  1(C) 

are  conformal  Since  every  link  of  P„  and  of  each  of  the  S  ’s  is  directed  in  the 
'  _  ‘  ‘  -  // 

.  same  sense  as  the  corresponding  link  of  P .  .  This  completes  the  proof. 


Although  this  result  is  needed  only  in  the  form  stated,  it  is  true  that  in 

general  m(P)  =  m(C)(+  Z’,  where  C  is  the  same  A/B  chain  as  above  and  Z’  is 

Hx  '  * ,  ■’  * 

obtained  by  multiplying  all  flows  of  Z  by  m. 


Lemma  1  is  used  to  establish  the  following  general  result  concerning  the 
decomposition  of  any  A/B  flow  pattern.  $■ 


Lemma  2: 

0 


If  F  is  an  A/B  flow  pattern  of  value  k  >  0, 
then 

F  =  (+)  1(C 2)(+) .  .  .01(Ck)0Z, 


(0 


where  the  C^s  are  appropriate  A/B  chains, 
Z  is  a  zero  flow  pattern,  and  Z  and  the 
l(C.)’s  are  conformal. 


o 


li 


i* 


The  proof  involves  a  constructive  procedure  that  produces  the  desired 


decomposition  and  makes  use  of  the  evident  fact  that  if  Y(M)  >  0  for  a  node  M 


there  is  at  least  one  link  M,  N  having  positive  flow,  i.  e. ,  there  is  a  link  having 
flow  directed  away  from  M.  This  property  is  true  for  any  assignment  of  flows 

•  O  ..  >> 

to  the  links  of  a  network,  whether  the  assignment  constitutes  a  flow  pattern  in 
the  sense  of  this  paper  or  not.  (In  the  course  of  the  proof,  F  undergoes  a  series 
of  modifications  that  at  intermediate  stages  determine  flow  assignments  that  are 
not  patterns. ) 


Proof:  Since  Y(A)  =  k  >  0,  a  link  A,  N  having  positive  flow  exists.  Let 


F^  denote  the  system  of  flows  obtained  from  F  by  reducing  the  flow  in  A,  by 


one  unit  (and  consequently  increasing  the  flow  in  N. ,  A  by  one  unit)  and  leaving 
all  other  link  flows  unchanged.  Then  Y(A)  =  k  -  1  >  0,  and  Y(N1)  =  1,  relative 


to  F^.  Hence  there  is  a  link  f^Ng  with  positive  flow.  Let  Fgbe  the  system  of 


flows  obtained  from  F1  by  reducing  the  flow  in  N^Ng  by  one  unit.  Then  Y(N^).  =  0 


again,  but  now  Y(Ng)  -  1  relative  to  Fg  (unless  Ng  is  B)v  So  long  as  each  suc- 


» 


cessive  Nj  is  not  B,  another  link  with  positive  flow  can  be  found,  be¬ 


cause  either  is  an  intermediate  node,  in  which  case  Y(N^)  =  1  relative  to  F^ , 


or  else  N.  =  A.  In  the  latter  case  Y(N.)  =  k.  (This  is  true  because  the  preceding 

■'])  3  3 


link  was,  of  the  form  N^,  A,  and  when  a  unit  of  flow  Was  subtracted  from  this 


link  in  passing  from  F^  to-5y  Y(A)  was  transformed  back  from  k  -  1  to  k. ) 

In  any  case  Y(N.)  >  0  if  N.  4=  B.  However,  this  procedure  must  terminate  since 


(i  3  ,v  3  "  *  .  <?- 

at  each  stage  a  unit  of  flow  is  subtracted  from  some  link  having  positive  flow. 


Thus  N,  =  B  at  some  stage  J.  But  then 


>  Tl 


Cc-  ,, 


F  =  FJ©1(A,N1,  .  .  .  ,  NJ_1,  B) 
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since  the  last  term  simply  represents  the  collection  of  unit  link  flows  removed 
from  F  to  obtain  F_.  Moreover  FT  is  an  A/B  flow  pattern  of  value  k  -  1,  being 
the  difference  between  F  and  a  unit  A/B  path  flow.  If  k  -  1  >  0  the  above  pro- 


cedure  can  be  applied  to  FT,  decomposing  it  into  an 

J 


flow  pattern  of  value 


// 


k  -  2  and  a  second  unit  A/B-  path  flow.  Ultimately 

;  F  =  Z’® 1^)©  1(P2)©  .  .  .  ©  l(Pk) 

where  Z  ’  is  a  zero  flow  pattern,  the  P.’s  are  A/B  paths,  and  the  nonzero  link 

0  *  ,  - 

flows  in  each  of  these  component  patterns  are  directed  the  same  way  as  in  F. 
Applying  lemma  1,  each  1(B.)  can  further  be  expressed  as  !(€.)(+,  Z.,  where 
C.  is  a  chain,  Z.  a  zero  flow  pattern,  and  all  link  flows  of  1(C.)  and  Z.  are 


directed  as  in  P. .  Thus 

o  1 


f  =  z©i(c1)©1(C2)©.  .  .'±:i(ck)’ 

J 

■  \ 

where  Z  =  ZT©  Z'  ©-  Z2©  .  .  .  ©  Zk  and  the  component  flow  patterns  are  all 

o 

conformal ,  This  completes  the  proof. 


1 


Feasible  Flow  Patterns  o 

A  flow  pattern  in  a  capacitated  network  is  said  to  be  a  ’’feasible  flow 

o 

pattern”  if  the  flow  in  every  link  is  feasible,  i.e. ,  if 


~c(N,  M)  <  f(M,  N)  <  c(M,  N) 


for  every  directed  link  M,  N.  There  are  only  a  finite  number  of  distinct  feasible 

\'v  •  * 

flow  patterns  that  can  be  associated  with  a  given  capacitated  network.  Letting 

V  . .  ‘ 

Vi  '  .  ,-v  c 

C  =  max  c(M,  N)  and  L  be  the  number  of  undirected  links,  there  are  at  most 


o 


2C  +  1  distinct  flows  possible  in  any  given  link,  viz.  ,  flows  of  1,  2,  .  .  .  ,  G 

L 

units  in  either  direction  or  zero  flow.  Thus  there  are  at  most  (2C  +  1)  ways 
to  assign  flows  to  all  lirks.  ♦  Of  course  most  of  these  assignments  are  not  flow 
patterns  in  the  sense  of  this  paper. 

-  *  '  c 

’  .  ~Z  (s' 

Lemma  3,  which  establishes  a  useful  result  concerning  feasible  flow 
patterns  states:  _ 


Lemma  3: 


If  F  is  a  feasible  A/B  flow  pattern  and 
Gl’  G2’  *  *  *  *  G.  are  conformal  A/B 
flow  patterns  such  that  F(+  G^©G2(+: . .  Q 
is  feasible,  then  F(+,  S  is  also  a  feasible  " 
A/B  flow  pattern  if  S  is  the  sum  of  any 


subset  of  the  G^’s. 


9  •  •  •  f 


Proof:  Let  f,  g1,  g^,  ...  ,  and  s  denote  the  corresponding  flows  in 
an  arbitrary  directed  link  M,  N.  Since  the  g^’s  are  all  nonnegative  or  all  non¬ 
positive,  one  of  the  following  relations  is  true: 

f<f+s<f+gi  +  g2+.  .  .  +  g,  '  ^  , 


gl  +  g2  +  *  •  *  +  gk -f+  8  -f 


♦More  precisely  there  are  exactly  Hf  c(M,  N)  +  c(N,  M)  +  1]  ways  to 
assign  flows  to  all  links  of  a  network,  where  the  product  ranges  over  all 
undirected  links  (M,  N)  . 


The  feasibility  assumptions  of,  the  lemma  also  imply  that  both  the  following 

© 

relations  hold: 


and 


-c(N,  M)  <  f  <  c(M,  N) 


€ 


-c(N, M)  <f  +  gx  +  g2  +  .  .  .  +  gk<  c(M, N) 

:  ",  -  JJ 


It  follows  that  ~c(N,  M)  <  f  +  s  <  c(M,  N);  I,  e. ,  the  flow  in  M,  N  corresponding  to 


the  flow  pattern  F(+  S  is  feasible.  But  since  M,  N  was  arbitrary  it  follows  that 
this  flow  pattern  is  feasible.  This  completes  the  proof. 


Using  lemmas  2  and  3  the  self-evident  fact  can  be  rigorously  established 


Jh 


that  if  a  feasible  A/B  flow  pattern  of  value  k  exists  in 'a  capacitated  network, 

/  \ 

feasible  A/B  flow  patterns  having  values  0,  1,  2,  ...  ,  k  exist.  For  if  Fk  is 
feasible  and  V(Fk)  =  k,  F^  can  be  expressed  as  Z(+  1(0^  (£  .  *  .  0  l(ck)  ac¬ 
cording  to  lemma  2.  If  Fj  is  defined  as  Z0 1(0^0.  .  .01(0.)  for 

f. 

i  —  1,  2,  .  .  .  ,  k,  and  FQ  =  Z,  according  to  lemma  3  F.  is  feasible  for 
1  =  0,  1,  .  .  .  ,  k.  Moreover,  since  Z  and  the  l(C.)’s  are  conformal,  the  F.’s 
are  nested  flow  patterns  in  the  sense  that  the  magnitude  of  flow  in  each  undi¬ 
rected  link  is  monotonic  nondecreasing  and  its  direction  remains  the  same,  as 
i  increases  from  0  to  k.  •'  „  0  ^  - 

A  "maximal  A/B  flow  pattern"  is  a  feasible  A/B  flow  pattern  F  whose 

(T^  °  ,  VX 

value  Y(F)  is  a  maximum.  Since  there  are  only  a  finite  number  of  feasible  flow 
patterns  possible  in  a  capacitated  network,  a  maximal  A/B  flow  pattern  exists 
for  every  pair  of  distinct  nodes  A  and  B.  (Note  that  the  flow  pattern  that  consists 
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o 


of  zero  flow  in  every  link  is  necessarily  feasible  and  can  be  considered  an  A/B 
flow  pattern  for  any  distinct  nodes  A  and  B. )  The  value  of  a  maximal  A/B  flow 
pattern  is  called  the  "A/B  capacity  of  the  network. " 


l 


0 


"  Cost"  of  a  Flow  Pattern 

c  Consider  next  an  A/B  flow  pattern  F  in  a  weighted  network.  F  is  associated 
with  the  nonnegative  integer  T(F)  defined  by  the  relation: 

T(F)  =  J||f(M^)|  «-  P  (M,  N)] 

S 

Where  S  is  an  arbitrary  enumerating  set  and  l’(M,  N)  is  defined  as 

1»(M,  N)  =  1(M,  N)  if  f(M,  N)  >  0  -** 


and 

1»(M,  N)  =  1(N,M)  if  f(M,N)  <  0 

/ 

Thus  T(F)  is  obtained  by  multiplying  the  flow  in  each  link  by  the  length, 

// 

measured  in  the  direction  of  the  flow,  and  summing  over  all  links  of  the  net- 

l  1 


work.  If  1(M,  N)  is  interpreted  as  the  cost  of  moving  one  unit  of  flow  from  M  to 
_  \ 

N  through  M,  N,  T(F)  is  simply  the  total  cost  associated  with  the  flow  pattern  F. 

h 

;u 

T(F)  is  necessarily  nonnegative.  ^  . 

C  “  'j/  :  .  - 

The  following  lemma  expresses  an  important  relation  between  the  costs 

,  Xyj 

of  certain  flow  patterns. 


Lemma  4:;  If  F,  G,  and  H  are  A/B  flow  patterns  and  G  and 

H  are  conformal,  then  \  \ 

.  .  ‘  '  V. 

T(F0G0H)  f  T(F@G)  >  T(F@H)  -  T(F) 
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Proof:,  Let  Q  =  T(F($G©H)  -  T(F®  G)  -  T(F@H)  +  T(F).  It  is  desired 
to  prove  that  Q  >  0  under  the  assumptions  of  this  lemma.  Let  S  be  an  enumerating 

V  O  -  .  _ 

set  chosen  in  such  a  way  that  f(C,  D)_>  0  for  all  C,  DeS,  where  f(C,  D)  denotes  the 

flow  in  C,  D  corresponding  to  pattern  F.  In  the  following  discussion  an  arbitrary 

_  .1  ^  _ °  . 

C,  DeS  is  fixed.  For  conciseness  the  flows  in  C,  D  corresponding  to  the  three 

.  9  '  _  c>  .. 

patterns  are  denoted  by  f ,  g,  and  h  and  the  lengths  1(C,  D)  and  1(D,  C)  by  1  and  1* 
respectively.  If  Q(C,D)  is  defined  by  the  relation 

-  Q(C, D)  =  |f  +  g  +  h|  xx  - [f  +  g |  x2  -  | f  +  h |  x3  +  |f  Jx4 

where  each  x.  is  1  or  1’,  depending  on  whether  the  corresponding  expression 
within  absolute  value  signs  is  positive  or  negative,  Q(C,  D)  represents  the  con¬ 
tribution  to  Q  corresponding  to  C,  D  and  S  Q(C,  D)  =  Q.  Thus  it  is  sufficient 
(\  J  \  C,DeS 

to  show  that  each  Q(C,  D)  >  0. 

One  of  the  following  sets  of  relations  holds; 

f  >  0,  g  >  0,  h  >  0  (3) 


It  remains  to  investigate  the  ease  when  relation  4  holds  and  f  +  g  +  h  <  0. 


In  this  case  Q(C,  D)  reduces  to  one  of  the  following  expressions,  depending  on  the 
relative  magnitudes  of  f,  g,  and  h:  ^ 

<i> 

-(f  +  g  +  h)l'  -  (f  +  g)l  -  (f  +  h)l  +  f  •  1 

,  a 

-(f  +  g  +  h)F  +  (f +  g)F  -  (f  +  h)  1  +  f  •  1 
~(f  +  g  +  h)F  -  (f  +  g)l  +  (f  +  h)F  +  f  •  1 
-(f  +  g  +  h)l*  +  (f  +  g)l*  +  (f  +  h)F  +  f  •  1 


\\ 


..  A 


When  simplified  these  reduce  to  the  expressions 

-(f  +  g  +  h)  (F  +  1) 

O 

-h(F  +  1) 

*  -k(!‘  + 1) 

f(F  +  1) 

These  are  all  nonnegative  since  f  g  +  h  <  0,  F  +  1  >  0,  and  the  conditions  of 

relation  4  hold.  Thus  Q(C,  D)  >  0  in  all  cases.  This  completes  the  proof.]' 

o 


a 


Thus  the  increment  of  cost  resulting  from  adding  H  to  F®  G  is  at  least 
as  great  as  that  resulting  from  adding  H  to  F  alone.  (This  result  fails  to  hold 
if  G  and  H  are  not  conformal.) 


Ideal  Flow  Patterns  ^  .  •  ; 

O  '' 

Finally  consider  a  network  that  is  both  capacitated  and  weighted,  and 
let  A  and  B  be  two  fixed,  distinct  nodes  of  the  network.  If  a  feasible  A/B  flow 

pattern  of  value  k  exists]  then  one,  call  it  F,  exists  such  that  T(F)  is  minimum, 

^  o  ■■  °  .  (  ^  .  ° 

since  the  feasible  A/B  patterns  of  value  k  are  finite" in  number.  Such  a  pattern 


is  termed  an  "ideal  A/B  flow  pattern  of  value  k. "  If  the  network  is  interpreted 
as  a  mathematical  model  of  a  transportation  network,  an  ideal  A/B  flow  pattern 
of  value  k  can  be  interpreted  as  a  pattern  that  moves  k  units  of  flow  per  unit 
time  from  A  to  B  at  minimum  cost.  Thus  ideal  flow  patterns  are  in  this  sense 
"best"  flow  patterns.  °  / 

The  following  lemma  establishes  an  importation  relation  between  certain 

o 

ideal  flow  patterns. 


Lemma  5: 


If  Fj  y  F.,  and  Fi+1  are  ideal  flow  patterns 
whose  values  are  i-1,  i,  and  i+1  respectively, 


T(Fi+1)  -  T(F.)_><T(F.)  -  T(F._1) 


Proof:  F.+1©F._1  is  an  A/B  flow  pattern  of  value  2.  As  such  it  can 
(according  to  lemma  2)  be  expressed  as  1(C1)@1(C2)@  Z  where  and  are 
A/B  chains;  Z  is  a  zero  flow  pattern;  and  1(0^),  1(C2),  and  Z  are  conformal. 
Since  Fj_.^  and  !  F.^1@1(C1)©1(C2)@  Z1  are  feasible  patterns  (the  latter  being 
simply  Fi+1),  it  follows  (according  to  lemma  4)  that  T(F.+1)  >  T(  F.^©  lfC.^) 
V  ©^(Cg)]  +  T(F.-1®  Z)  -  T(F.-1).  But  F.__1©  Z  is  an  A/B  flow  pattern  of  value 
i-1  and  is  feasible  (according  to  lemma  3).  Hence  T(F.  ^®  Z)  -  T(Fi_1)  >  0 
since  is  ideal.  It  follows  that 


Applying  lemma  4  again 


T(F.+1)  >  T|  Fi>_1®l(C1 


>  n 'Fi_1©l(C1)  +  T(F._1©1(G2)]  -  TfF^j) 


But  1(C|)  arid  F^_^v+  1  (C2)  are  A/B  flow  patterns  of  value  1  and  are 

feasible  (according  to  lemma  3).  Since  F.  is  ideal,  T[  Fi_^01(C1)]  >  T(F.) 
and  T[  Y.^&nC^  >  T(F.).  Hence 

T(Fi+1)  >  2T(F.)  -  T(F._1)  \ 


But  this  is  equivalent  to  the  inequality  that  Was  to  be  established. 


A  sequence  of  real  numbers  aQ,  a^,  a2 


9  •  •  -•y,  r 


*k 


is  said  to  be  convex 


if 


w  .4 


a..  -  -  a.  >  a.  -  a.  , 
i+I  i  -  i  i-l 


for  i  =  1,  2,  ...  ,  k-1.  If  Fq,  F^,  Fg^v  •  »  F^  are  ideal  A/B  flow  patterns 
of  values  0,  1,  2,  .  ,  k  respectively,  where  k  is  the  A/B  capacity  of  the  net¬ 

work,  it  follows  that  the  sequence  T(Fq),  T(F1),  . .  ,  T(Ffc)  is  a  convex 
sequence  of  nonnegative  real  numbers,  the  first  of  which  is  zero.  I  The  flow  ' 


pattern  Fft,  which  assigns  zero  flow  to  every  link,  is  clearly  ideal,  and 

V  0  '  °  \$l 


T(F0)  =  0.] 

In  a  certain  sense  the  convex  sequence  T(Fq),  T(Fj),  .  .  .  ,  T(F^)  can 
be  considered  as  the  "A/B  cost  profile"  for  a  given  network.  The  remainder  of 
this  paper  is  devoted  to  developing  a  general  method  for  determining  the  A/B  cost 


profile  for  a  given  network  and  given  nodes  A  and  B,  and  for  finding  specific 


V 


ideal  flow  patterns  possessing  these  minimum  costs  .  The  problem  is  precisely 
defined  in  the  next  section. 

As  a  matter  of  interest,  note  that  given  any  finite  convex  sequence  of 

*  0-  " 

nonnegative  real  numbers  aQ,  a^,  a2,  . 


,  an  a  network  exists  having  this 


23 


sequence  as  its  A/B  cost  profile.  In  fact  the  network  of  Fig.  2  has  the  desired 
properties.  In  this  figure  the  numbers  indicate  the  lengths,  in  both  directions, 
of  the  corresponding  links.  All  link  Capacities  are  assumed  to  be  one  unit.  If 
1  <  i  <  n,  the  (in  this  case  unique)  ideal  A/B  flow  pattern  of  value  i  consists  of 
assigning  one  unit  of  flow  from  A  to  B  in  each  of  the  paths  A,  N.,B,  j  =  1,  2,  . . . ,  i, 

--  V,  ^  ^ 

and  assigning  zero  flow  in  the  remaining  links.  Clearly  T(F)  =  for  this  pattern. 

(The  intermediate  nodes  N. ,  N0,  .  .  ~ N  are  introduced  merely  to  conform 

Jl  a  n 

•  -  ° 

with  the  requirement  that  at  most  one  link  join  two  distinct  nodes. ) 

G  '  ij 

//  *  . 


o 


Fig.  2 — A  Network  Having  A  Given  A/B  Cost  Profile 


o 


STATEMENT  OF  PROBLEM  AND  SOLUTION  IN  PRINCIPLE 


formal  Statement  of  problem 

In  the  remainder  of  this  paper  the  network  under  consideration  is 

assumed  to  be  connected,  weighted,  and  capacitated  and  to  have  the  property 

// 

that  for  every  cycle  S,  1(S)  >  0. 

The  problem  of  immediate  concern  can  be  precisely  stated: 

Given  a  network  of  the  above  type,  and 
two  distinguished  nodes  A  and  B,  produce  a 
procedure  that  will 

i 

Qa)  determine  the  A/B  capacity  k  of 

Ji 

■  t  :  : 

the  network; 

^  / 

(b)  determine  a  specific  ideal  A/B  flow 

pattern  F.  of  value  i,  for  i  =  1,  2,  .  ,  k;  and 

e>  1 


t>  '  ' 


0  (c)  always  terminate  in  a  finite  number 

of  steps. 

.  ?/  Network  flow  problems  of  various  types  have  been  the  subject  of 

considerable  study  and  a  number  of  procedures  having  property  (c),  that  can 

C 

be  used  to  solve  problem  (a),  or  (b)  for  any  specified  value  of  the  flow  pattern , 
exist  in  the  literature.  Reference  is  made  to  some  of  these  after  the  present 
procedure  is  defined  and  its  validity  established.  sO 

In  this  section  a  procedure  solving  the  above  problem  in  principle  is 
given,  accompanied  by  the  proofs  required  to  establish  its  validity.  The 


25 


If 


'*>x 


development  is  incomplete  in  one  respect,  however.  The  assumption  is  made 
that  it  is  possible  to  design  an  algorithm  for  finding  A/B  chains  of  a  certain 
type.  The  next  section  of  the  paper  produces  this  algorithm,  which  completes 

C 

the  development. 


AN  EXAMPLE 


The  following  example  motivates  the  central  concept  involved  in  the 
procedure,  viz. ,  the  effective  length  of  a  chain. 

Consider  the  network  shown  in  Fig.  3a.  Here  A,  N2>  N_,  N^,  and  B 
represent  the  nodes,  and  there  are  seven  links  as  indicated.  The  numbers  above 
and  below  each  link  represent  its  length  and  capacity  respectively.  (These  apply 
to  both  directed  links:  i.e. ,  the  network  is  fully  symmetric. )  An  ideal  A/B 

O  •  '  -  ‘  ' 

flow  pattern  of  value  0  is  of  course  the  pattern  FQ,  which  assigns  zero  flow  to 
every  link.  The  unit  A/B  chain  flow  F^  =  1(A,  N2,  Ng,  N^,  B)  is  an  ideal  pat- 

O 

tern  of  value  1.  It  happens  to  be  the  only  one  for  this  network.  Note  that 
T(Fj)  =  4  and  that  saturates  the  directed  links  A,  Ng,  N9,  Ng,  and  Ng,  N^. 

One  possible  way  to  augment  F^  and  obtain  an  A/B  pattern  oT value  2t  G2» 

ty,  _  ’’ ! 

would  be  to  add  1(A,N^,  B)  to  F^.  The  resulting  pattern  is  shown  in  Fig.  3b, 
where  the  numbers  and  arrows  indicate  the  magnitudes  and  directions  of  link 
flows.  But  G2  is  not  ideal  since  T(G2)  -  11,  whereas  T(F2)  =  10,  where  F2  is 
the  feasible  A/B  flow  pattern  shown  in  Fig.  3c.  Fg  is  ideal,  and  there  is  no 
other  ideal  pattern  of  value  2  in  this  network.  F2  is  also  maximal  since  clearly 

O 

link  N^,  B  limits  flow  to  two  units  at  most.  Thus  the  cost  profile  for  this 
network  is  0,  4,  and  10. 


o 


Now  considering  the  statement  of  the  problem,  especially  property  .(b),, 
it  is  highly  desirable  that  the  procedure  developed  involves  using  F.,  once  it  is 

found,  as  the  starting  point  for  constructing  Fi+1  rather  than  constructing  F.+1 

})  . 

independently.  1  &- 

If  Fg  is  to  be  obtained  by  augmenting  F^^  in  the  example,  the  pattern  that 
must  be  added  to  F.^  is  F g£  F^  =  1(A,  Ng,  Ng,  N^,  B) .  It  is  not  coincidental  that 
Fg@  Fi  is  a  unit  A/B  chain  flow.  It  is  shown  below  that  if  Fi  is  ideal  and  non- 
maximal,  an  ideal  pattern  F.+1  can  always  be  obtained  by  adding  an  appropriate 
unit  A/B  chain  flow.  There  may  be  other  ways  also,  of  course.  In  the  case  of 
the  present  network,  since  Fj  and  Fg  are  unique,  it  follows  that  Fg0  F^  had  to 
be  a  unit  A/B  chain  flow. 

Once  the  above  result  is  rigorously  established,  viz. ,  that  it  is  sufficient 
to  consider  patterns  of  the  form  F.  01(C)  in  generating  F.+1,  it  remains  to  es¬ 
tablish  a  simple  criterion  for  determining  which  chain  to  use  and  to  develop  a 
systematic  method  for  finding  a  chain  satisfying  the  criterion. 

Consider  again  the  chains  involved  in  generating  Gg  and  Fg  from  F^  in 
the  example.  These  are  Cj  =  A,  N^,  B  and  Cg  =  A,  Ng,  Ng,  N^,  B  respectively. 

*  In  terms  of  length,  C1  is  preferable  since  1  (Cj)  -  7  whereas  l  (Cg)  =  8.  In 
terms  of  increments  of  T,  however,  T(Gg)  -  T(F^)  =  7  and  T(Fg)  w  T(Fj^)  -  6. 
Thus  Cg  is  effectively  shorter.  This  notion  is  now  made  general  and  precise. 


EFFECTIVE  LENGTH 


If  F  is  a  feasible  A/B  flow  pattern  with  associated  link  flows  denoted  by 


'zr? 


ij  f(M,  N),  where  M,  N  is  an  arbitrary  directed  link,  the  "effective  length  of 


M,  N  relative  to  F1 ,.  denoted  by  e(M,  N;F),  is  defined  as:  / 

Jj  ' 

^  '  f=  1(M7N)  if  o  <  f(M^N)  <  c(M7N) 

t  •. 

e(M^N;F)<  =  -1(Xm)  if  f(M7N)  <  0 
0  [=  «*  if  f(MTN)  =  c(M^N) 

If  C  denotes  any  chain  or  cycle,  the  "effective  length  of  C  relative  to  F’, 
denoted  by  e(C;F),  is  defined  as: 


e(C;F)  =  /  e(M,  N;F) 

e 

V 

Note  that  e(M,  N;  F)  =«  if  and  only  if  F  saturates  M,N,  and  consequently 
e(C;F)  <»  if  and  only  if  no  link  of  C  is  saturated.  This  leads  to  the  following 

O  ^ 

result: 


Lemma  6: 


If  F  is  a  feasible  A/B  flow  pattern  and  C  ^ 
is  an  A/B  chain  or  a  cycle,  F(+)  1(C)  is 


~  feasible  if  and  only  if  e(C;F)  is  finite.  ^ 

Proof:  A  link  M,  N  is  saturated  if  and  only  if  e(M,  N;F)  =°°,  Thus  if 


e(C;F)  is  finite,  no  link  of  C  is  saturated,  and  so  the  flow  in  each  link  of  C  can 

f\  '  , 

,  %  ° 

be  increased  by  at  least  one  unit.  Hence  F@1(C)  is  feasible  if  e(C;F)  is  finite. 
Conversely  if  F(+,1(C)  is  feasible  no  link  of  C  is  saturated,  and  e(C;F)  is  finite 
since  e(M,  N;F)  is  finite  for  M,  NeC,  This  completes  the  proof.  ^ 


♦For  convenience  "•*>"  is  adjoined  to  the  real  number  system.  The  only 
properties  of  this  element  that  are  required  are  (a)  if  R  is  any  real  number, 

R  <  «o;  (b)  it  is  not  true  that  »  <  «>;  (c)  R  +  «*  =  •;  (d)  «  +  «*  =  «;  and 
(e)  min{  •,  «°,  .  ,  ...  ,  «*}  =  •.  0  ^ 

’  ‘  .  '  "  v  /  . V;  '■  /"• 
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Lemma  7  is  a  formal  statement  of  the  fact  that  the  definition  of  effective 
length  satisfies  the  property  suggesting  its  introduction. 

Lemma  7:  If  F  is  a  feasible  A/B  flow  pattern  and  C  is 

an  A/B  chain  or  a  cycle  such  that  e(C;F)  is 
*  finite,  then  T[  F@  1(C)  ]  -  T(F)  =  e(C;F) . 

Proof:  Consider  any  M,  NeC.  If  f(M,  N)  >  0,  where  f  denotes  flow 
corresponding  to  F,  then  the  flow  in  M,  N  corresponding  to  F®  1(C)  is  f(M,  N)  +1 
and  T  is  increased  in  this  link  by  the  amount  1(M,  N)  .  If  f(M,  N)  <  0,  the  flow  in 
mTn  corresponding  to  F®  1(C)  is  f(M,  N)  +  1  <  0  and  that  in  N,  M  is  f(N,  M)  -  1  >  0 
thus  t  is  decreased  in  this  link  by  1(N,  M)  units  or  increased  by  -1(N,  M).  In 
either  case  the  increase  in  T  with  respect  to  M,N  equals  e(M,  N;F)  as  defined 
above.  Since  M,  N  was  arbitrary,;  the  result  follows . 

The  following  result,  which  is  an  immediate  consequence  of  lemmas  6 
and  7,  plays  an  important  part  in  the  next  section.  _,<<* 

<N 

fS 

Lemma  8:  >.  If  F  is  an  ideal  A/B  flow  pattern,  then 

O  ■-  ' 

e(S;  F)  >  0  for  every  cycle  SJ\  " 

**  yC  U 

Proof:  According  to  lemma  7  T[  F®1(S)]  '=  T(F)  +  e(S;F).  If  e(S;F)  <  0 
it  follows  that  Tt  F®  1(S)]  <  T(F).  But  F®  1(S)  is  a  feasible  A/B  flow  pattern 

according  to  lemma  6  and  has  value  V(F),  which  contradicts  the  assumption  that 

„  '  *  '  *  /'  •> 

F  is  ideal.0  t( 


CENTRAL  THEOREMS  \\ 

•>  ^  •  0 

The  set  of  maximal  A/B  flow  patterns  can  be  characterized  in  terms  of 

the  notion  of  effective  length  in  the  following  way. 

o  ^ 

Theorem  1:  A  feasible  A/B  flow  pattern  F  is  maximal  if 

and  only  if  e(C;F)  =  «  for  every  A/B  chain  C. 

*  _  ’  ~  v 

Proof:  If  e(C;F)  is  finite  for  some  A/B  chain  C,  then  F®  1(C)  is  feasible 
by  lemma  6,  and  thus  F  is  not  maximal.  Hence  if  F  is  maximal,  e(C;  F)  =  »  for 
all  A/B  chains.  Conversely  if  F  is  not  maximal,  a  feasible  A/B  flow  pattern  G 
of  value  V(F)  +  1  exists..  According  to  lemma  2  G®F  can  be  written  as  1(C)®  Z 
where  C  is  an  A/B  chain,  Z  is  a  zero  flow  pattern,  and  1(C)  and  Z  are  conformal. 
Since  F  and  F®  1(C)®  Z  =  G  are  feasible  patterns,  it  follows  from  lemma  3  that 
F®1(C)  is  feasible  also.  Thus  e(C;F)  is  finite  according  to  lemma  6.  Hence  if 
-  e(C;F)  =  00  for  all  A/B  chains,  F  is  maximal.  This  completes  the  proof. 

Since  there  are  only  a  finite  number  of  distinct  A/B  chains  in  a  network 
and  each  has  a  well-defined  effective  length  relative  to  any  given  feasible  A/B 
flow  pattern  F,  a  chain  C.  exists  that  minimizes  e(C;F) .  This  chain  is  referred 
to  as  a  "chain  of  shortest  effective  length  relative  to  F’ ,  or  as  an  "effectively 
shortest"  chain.  It  is  a  consequence  of  theorem  1  that  e(CjjF)  <  00  if  and  only 
if  F  is  nonmaximal.  ' 

c  ^  o 

An  ideal  A/B  flow  pattern  that  ii  nonmaximal  can  always  be  extended  to 
an  ideal  pattern  of  value  one  unit  greater  in  a  m  anner  described  in  the  next 
theorem. 


ITT 


Theorem  2;  If  F  is  an  ideal  nonmaximal  A/B  flow 

pattern  of  value  k  and  C  is  any  A/B  chain 
of  shortest  effective  length  relative  to  F, 

t , 

then  F@1(C)  is  an  ideal  A/B  flow  pattern 

W  '  ’  ° 

of  value  k  +  1. 

Proof:  Let  G  be  an  ideal  flow  pattern  of  value  k  +  1.  Then  according  to 

"  •  '  {'■  ■*'<>  ; 

lemma  2  gQ F  can  be  ejqjressed  as  1(G,)@  Z,  where  C'  is  an  A/B  chain,  Z  is  a 
zero  flow  pattern,  and  l(C')  and  Z  are  Conformal.  Since  F  and  F(+/1(C')®Z  are 
feasible,  lemma  4  asserts  that  the  following  inequality  holds:  // 

ft  '! 

T(G)  -  T[  Fg)l(C)]  >  T(F©Z)  -  T(F) 

But  F@Z  is  feasible,  according  to  lemma  3,  so  that  T(F@  Z)  -  T(F)  >  0  since 
F  is  ideal.  Thus  T(G)  >TI  F©1  (C1)].  Now  F@1(C)  is  also  feasible,  according 
tj  to  lemma  3,  and  has  value  k  +  1.  It  follows  that  F01(C)  is  also  ideal.  C  must 
minimize  e(C;F)  because  if  e(C”;F)  <  e(C’;F)  it  would  follow  from  lemma  7  that 
T[  F@1(CM)J  <  Tl  F@i(C’|  i  It  is  also  clear  from  lemma  7  that  if  C  is  any 
other  chain  with  e(C;F)  =  e(C’;F),  then  F©1(C)  is  an  ideal  A/B  flow  pattern 
whose  value  is  k  +  1 .  This  completes  the  proof. 


If  the  existence  of  a  terminating  algorithm  that  will  always  find  a  chain  of 
shortest  effective  length  relative  to  an  ideal  flow  pattern  F,  when  F  is  non- 

•  C  .  ‘  ii  v.  •  in 

maximal, and  will  give  a  suitable  indication  when  F  is  maximal  is  assumed, 
then  theorems  1  and  2  could  be  used  to  form  the  basis  for  a  procedure  for 
solving  the  stated  problem.  However,  this  Would  require  proceeding  in 
increments  of  one  unit  of  flow— a  severe  limitation  in  any  practical  problem.  ... 


o 


c 


c  l*j  ;  •  ■ 

Ns:  •  C 

C- 

Fortunately  the  convexity  property  established  in  lemma  5  can  be  used  to 
formulate  a  general  method  for  accelerating  the  process. 

Let  F,  be  an  ideal  A/B  flow  pattern  of  value  i  in  a  network  whose  A/B 
...  capacity  is  at  least  i  +  2.  Then  if  the  A/B  chain  CL  minimizes 

1  y 

e(C;F),  Fi+1  =  F.®  l(C^)  is  ideal  according  to  theorem  2.  Similarly 
Fi+2  =  F.+1@1(C2)  is  ideal  if  C 2  minimizes  e(C;Fi+1):  Applying  lemma  7 
the  following  equations  hold: 

\\  :  i  : 

5  T(Fi+2)  -  T(Fi+1)  =  e(C2;Fi+1) 

and 

•  0 

T(F.+1)  -  T(F.)  =  e(Ci;F.) 

(> 

But  lemma  5  asserts  that 


Now  suppose  that  e(CL ;  F. , . )  -  e(C,;F,).  Since  C«  minimizes  e(C;F.  ),  it 

1  1+1  1  1  ’  *  l+l 

follows  that  is  an  effectively  shortest  chain  relative  to  Fi+1  also,  so  that 
>-G1+2  =  Fi+101(C^)  =  Fi©2(C1)  is  an  ideal  A/B  pattern  of  value  i  +  2.  In 
general  once  a  chain  CL  that  minimizes  e(C;F.)  has  been  found  ideal  patterns 
Fj+1,  Fi+2,  .  .  .  can  be  generated  by  adding  1(0^  to  each  preceding  pattern. 
As  long  as  e(C1;Fi+k_1)  =  e(C1;F.),  it  will  be  true  that  Ft+k  =  F1+k_£®  l(Ci) 
=  Fj®k(C  j)  is  an  ideal  A/B  pattern  of  value  i  +  k.  Theorem  3  establishes  a 


0 


o 


criterion  for  determining  the  maximum  number  of  times  that  the  same  unit  chain 

£  .  '  '  ... 
flow  l(Cj)  can  be  added  to  an  ideal  flow  pattern  before  the  resulting  pattern 

ceases  to  be  ideal. 


J 


Theorem  3: 

'  )  ^ 


If  F  is  an  ideal  nonmaximal  A/B  flow  pattern 
and  C  is  any  A/B  chain  of  shortest  effective 

u 

length  relative  to  F,  then  F0j(C)  is  an  ideal 

A/B  flow  pattern  of  value  V(F)  +  j  ,  for 

j  =  1,  2,  .  .  ,  Q,  where  Q  =  min  q(M,  N;F) 

C 


and  q  (M,  N;F)  is  defined  as 


a(MN^F\  =  C<M*N1  ”  f<M’N>  if  W>N)  > 
q'JVA’  ■  ’ *  '  -f(M,  N)  if  f (Mr  N)  <  0 


) 


J  ’■ 


Proof:  From  theorem  1  it  is  known  that  e(M,  N;F)  is  finite  for  every 


M,  NeG.  Referring  to  the  definition  of  e(M,  N;F),  two  cases  must  be  considered. 


(a)  f (M,  N)  >  0.  In  this  case  e(M,  N;F)  =  1(M,  N).  In  fact  e[  M,  N;F0k(CJ 


V. 


=  1(M,  N)  so  long  as  f(M,  N)  +  k  <  c(M,  N);  ive. ,  so  long  as  k  <  c(M,  N)  -  f(M,  N). ' 


1(G)  can  be  added  to  F  c  (M  ,  N)  -  f(M,  N)  -  1  times  without  altering  the  effective 

T  ■■  ‘  ■  — ‘  O 

length  of  M,  N.  Thus  1(C)  can  be  added  to  F  one  more  time,  [i.e. ,  a  total  of 


q(M,  N;F)  time^  before  the  effective  length  changes.  In  this  case  M,  N  is 


saturated  and  el  M,  N;F@q(C)]  =  ». 

a  ' 


(b)  f(M,  N)  <  Oj}.  In  this  case  e(M,  N;F)  = -1(N,  M).  In  fact  el  M,  N;F@k(CJ 


=  -1(N,  M)  so  long  as  k  <  -f(M,  N).  Thus  the  effective  length  of  M,  N  changes  only 

V>‘  -  .  /  .  ,  °1.  :  ^  '  ' 

after  adding  1(C)  to  F  ,  -f(M,N)  times,  which  is  the  same  as  q(M,  N;F)  times. 
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In  this  case  the  flow  in  (M,  N)  is  then  reduced  to  zero,  and  the  effective  length 
of  M,  N  changes  from  ~1(N,  M)  to  1(M,  N). 

New  if  Q  is  taken  to  be  the  minimum  of  the  link  q’s,  it  follows  that  for 
j  =  1,  2,  .  .  Q-l  '  "J 

^  e|C;F®j(C)]  =  e(C;F) 

j!  , 

since  this  relation  is  true  for  every  directed  link  of  C.  But  it  follows  from  the 
remarks  immediately  preceding  the  statement  of  theorem  3  that  so  long  as  the 
effective  length  of  C  remains  constant  the  addition  of  1(C)  to  each  successive 

,s  1  .  ft 

pattern  generates  another  ideal  pattern.  This  completes  the  proof. 

PROCEDURE 

The  existence  of  an  algorithm,  which  is  referred  to  as  the  "effective1- 
length  algorithm"  ,  is  assumed.  As  input  it  requires  an  enumerating  set  S 
characterizing  the  network  configuration  and  an  ideal  A/B  flow  pattern  F.  of 
value  i,  described  by  Itemizing  the  flow  f.(M,  N)  for  all  M,  NeS.  As  output  it 
produces  the  effective  distance  Dj  from  A  to  B  relative  to  F..*  KDj  <  «>  (i.e. , 

N-  '  fr  0 

%  if  Fj  is  nonmaximal)  it  also  produces  a  specific  effectively  shortest  chain  C.. 

A  procedure  that  will  always  solve  the  stated  problem  can  now  be 

formally  presented.  The  procedure  involves  repeating  the  following  steps  as  ,, 

’  ’  !} 

(c).. 

♦"Effective  distance"  means  the  effective  length  of  an  effectively 
shortest  chain. 


o 


many  times  as  required  to  reach  a  maximal  flow  pattern.  At  the  start  take 


i  =  0  ahd  take  as  F^  the  ideal  A/B  flow  pattern  for  which  f^ (M,  N)  -  0  for  all 

'\s 

M^NeS.  c  a  " 

(a)  Apply >ithe  effective  distance  algorithm  to  the  network  with  flow 
pattern  F. .  If  D.  =•'*>•  (i.e. ,  if  F.  is  maximal),  then  terminate  the  procedure 
because  the  complete  solution  to  the  problem  is  obtained.  Otherwise,  perform 
step  b . 

c  ,  o 

(b)  Let  C;  be  the  effectively  shortest  chain  found  in  a.  Sean  the  links 

1)\  ■  _  ..  \ 

of  C.  and  determine  0.  ~~  min  J  q(M,  N);F.]  as  defined  in  theorem  3> 

O  '  1  ^ 

(c)  Construct  an  ideal  A/B  flow  pattern  F:+^  of  value  i  +  k  as  follows7' 

for  k  '=  1,  2,  .  .  .  ",  Q.: 

#  1  *  .  . 

For  every  M,  NeS  define  f.+^  as 

<T-  V  _ 

(1)  f.  +  kif  M.NcC., 

(2)  f.  -  k  if  N,  MeC. 

(3)  f  otherwise 
*  1 

Then  repeat  (a)  with  i  +  Q.  replacing  i  and  F.+q  replacing  F.. 

Theorem  3  states  that  the  flow  patterns  obtained  in  this  way  are  ideal, 

and  according  to  theorem  1  the  process  will  terminate  when  and  only  when, 

after  step  c  is  completed,  the  resulting  F.  Q  is  maximal. 

&  1  u 

<>C’  ,  v  . 

RELATION  TO  OTHER  PROCEDURES 


Maximizing  Flow.0  The  principle  of  obtaining  a  maximal  flow  pattern  by 

O  .. 

adding  a  succession  of  unsaturated  A/B  chains  is  well  known.  Dantzig  and 


c 


Fulkerson  present  an  efficient  hand-computing  scheme  using  this  principle, 
which  is  applicable  to  any  capacitated  network1  with  symmetric  capacities  [  except 
that  c(M,  N)  =  0  if  M  ~  B  or  N  =  A]  *•  Any  procedure  for  determining  unsaturated 
A/B  chains  could  be  made  the  basis  for  a  flow-maximizing  process.  The  present 
procedure  is  a  special  instance  of  this  approach  since  it  selects  from  all  A/B 
chains  C  one  that  minimizes  e(C;E). 

I 

Minimizing  Cost  for  Stated  Flow.  If  k  is  a  nonnegative  integer  that 

V  '  ;  '  ’  '  ,  ■  ' 

does  not  exceed  the  A/B  capacity  of  the  netwoil ,  the  problem  of  finding  a  feasi¬ 


ble  A/B  flow  pattern  F  such  that  V(F)  =  k  and  T(F)  is  minimum  is  a  special  case 
of  the  general  linear-programing  problem.  Thus,  in  theory,  a  technique  such 
as  the  simplex  method  could  be  employed.  However,  in  complex  networks  the 
number  of  variables  and  constraints  is  so  large  as  to  make  such  a  general 

approach  infeasible.  \s  f  • 

:  \  ...  '■  :  . 

This  problem  is  also  a  special  case  of  the  class  of  network-flow  problems 
known  as  ’’capacity-constrained  transshipment  problems .  ”  Fulkerson2  describes 
this  type  of  problem  and  shows  that  it  is  equivalent  to  an  appropriate  Hitchcock 
problem.  Thus  any  procedure  for  solving  the  Hitchcock  problem  could  also  be 
employed. 

■  ^  o 

Finding  a  Complete  Set  of  Ideal  Flows.  The  problem  of  finding  a  family 
of  ideal  A/B  flow  patterns,  one  for  every  feasible  value  of  flow,  is  related  to  a 

C> 

dynamic  problem  posed  and  solved  by  Ford  and  Fulkerson  of  the  RAND  Corpora- 
3 

tion.  That  problem  deals  with  maximizing  the  total  flow  arriving  at  a  node  Be 
by  the  end  of  T  time  periods  assuming  there  is  no  flow  in  trans'l;  at  time  T  =  0. 

-rJ'  •-  \\  O  .  -■ 

(The  lengths  of  links  are  assumed  to  be  positive  integers  and  are  interpreted 


as  the  number  of  time  periods  required  to  traverse  the  links.  )  The  procedure 

//'  '  • 

presented  finds  a  family  of  solutions — one  for  every  value  of  T.  (These  > 
solutions  remain  the  same  after  a  certain  value  of  T  is  reached. )  The  pro¬ 
cedure  involves  finding  static  flow  patterns  in  the  network,  and  it  is  shown  that 
these  flow  patterns  are  ideal,  at  least  when  T  is  sufficiently  large. 

Decomposing  Solutions  into  Routes .  Although  the  procedure  described 
in  this  paper  does  not  include  a  specific  method  for  expressing  the/ideal  flow 

(t  i 

\\ 

patterns  as  the  sum  of  conformal  chain  flows ,  *  it  is  recognized  that  this  is 

N\ 

'  '  o  •  V  <■'•■ 

desirable  if  one  wants  .an  operational  plan  for  routing  shipments  at  minimum 

''  {}  c> 

cost.  The  construction  employed  in  the  proof  of  lemma.  2  could  be  formalized 

into  an  algorithm  for  producing  such  a  composition,  t  However,  efficient 

schemes  for  decomposing  flows  exist.  Specifically  the  labeling  process 

3  : 

termed  "Routine  H"  by  Ford  and  Fulkerson  performs  this  function  efficiently. 


♦The  example  given  at  the  beginning  of  the  present  section  shows  that 
the  chain  flows  used  to  synthe  ttze  an  ideal  pattern  are  not  necessarily  conformal. 

t Referring  to  lemma  2  the  residual  flow  pattern  Z  is  necessarily  the 
pattern  that  is  identically  zero.  This  can  be  shown  to  be  a  consequence  of  the 
assumption  that  every  cycle  C  has  positive  length.  If  Z  were  not  identically 
zero,  the  pattern  F  being  decomposed  could  not  be  ideal. 


V.' 


EFFECTIVE-LENGTH  ALGORITHM 


STATEMENT  OF  THE  ALGORITHM 


Let  F  be  an  ideal  A/B  flow  pattern  in  a  network,  with  link  flows  denoted 

i  by  f(M,N).  Consider  the  links  of  the  network  to  be  arranged  in  a  sequence 

-  (M.,N.),  i  =  I,  2,  .  .  .  ,  L.  ,The  algorithm  consists  of  the  assignment  and 

revision  of  certain  quantities,  or  "labels",  associated  with  the  nodes.  It  is 

patterned  after  the  labeling  procedures  employed  by  Ford  and  Fulkerson  for 

4 

solving  various  network-flow  problems.  Ford  describes  a  labeling  process 

concerned  with  finding  an  A/B  Chain  C  that  minimizes  1(C)  =  2 1(M,N).  The 
O  C 

present  algorithm  is  an  adaptation  of  that  one,  where  the  necessary  modifications 

were  made  to  minimize  2  e(M,N;F)  rather  than  2 1(M,N).  In  addition  infor- 

C  df  , 

/  mation  is  carried  along  that  enables  one  to  readily  determine  a  specific  chain 

that  minimizes  2  e(M,N;F).  This  is  achieved  by  recording  certain  approach 

C  _  .  1 

o  0 

links  in  the  same  manner  as  the  labeling  process  of  Ford  and  Fulkerpon, 

/  ‘  '  ‘  \ 

Formal  description  of  the  effective-length  algorithm  follows. 

}  " 

(a)  Associate  with  node  A  the  label  D(A)  =  0,  and  let  D(N)  «  for  all 
v other  nodes. 


(b)  Considering  each  link  (M.,N.)  in  turn,  perform  the  following 

■  1  1  /! 


adjustments: 


o  (1)  If  D(N.)  >  D(M.)  +  e(Mf,N.;F),  replace  D(N.)  by  this 

1  ,  1  11  :  VV.r>  1 

smaller  quantity  and  record  M.,N^  as  the  approach  ih.k  currently  associated 


ith  node  N. . 
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(2)  If  D^Mj)  >  D(N.)  +  e(N.,M.;F),  replace  D(M.)  by  this  smaller 
quantity  and  record  N. ,  M;  as  the  approach  link  associated  with  node  M. ... 

(3)  If  neither  of  the  foregoing  inequalities  holds,  make  no  change. 

f'.  v  \  ... 

Repeat  step  b  until  the  stage  is  reached  when  3  applies  for  every  link 
(ML,  N. )  of  the  network.  V.,  !  ° 

O  .  o  ); 

It  is  shown  below  that  this  stage  is  necessarily  reached  in  a  finite 

•%  •  . 

number  of  steps,  that  at  this  time  D(B)  =?  min  e(C;F)  for  all  A/B  chains  C,  and 

,  i  c  . 

that  the  final  set  of  approach  links  can  be  used  to  determine  a  specific  effectively 
shortest  A/B  chain  C,  if  min  e(C;F)  <  00  ,  i.  e. ,  if  F  is  nonmaximal. 

DERIVATION  OF  PROPERTIES  OF  THE  ALGORITHM 

(M 

The  next  five  results  (lemmas  9  to  13)  state  certain  relations  that  exist 
between  finitely  labeled  nodes  and  approach  links  at  any  stage  of  the  algorithm  . 

i ! 

Assume  that  F  is  an  ideal  A/B  flow  pattern  in  a  network  and  that  after  starting 
to  apply  the  rules  of  the  algorithm  the  process  is  interrupted  at  an  arbitrary 
time.  The  labels  and  approach  links  described  in  the  statements  and  proofs  of 
these  lemmas  are  assumed  to  be  those  that  exist  at  the  time  of  interruption 
except  where  otherwise  noted. 


)>  - — 

Lemma  9:  IfM,N  is  an ‘approach  link,  the  following 

relation  holds:^ 


c-j- 


V' 


Proof:  Let  D(M)  denote  the  label  of  node  M  at  the  time  it  was  used  to 


assign  the  present  label  to  node  N.  Then 


D(N)  =  D(M)  +  e(M,  N;F) 


as  a  result  of  the  relabeling  rules .  But 


since  labels  are  never 


increased  by  the  process.  The  inequality  of  the  lemma  follows. 


Lemma  10: 


The  set  of  approach  links  do  not  contain  a 
subset  that  is  a  loop. 


Proof:  Suppose  Np  Ng,  .  .  .  ,  (N^  =  is  adbop,  where  each 

N. ,  N.+1  is  an  approach  link.  This  loop  is  necessarily  a  cycle.  (Since  a  node 

never  has  two  approach  links  associated  with  it,  the  terminal  nodes  of  the  links 

//  ■  ■ 

are  all  distinct. )  According  to  lemma  9  then 


+  e(Ni.r  Nt;F) 


for  i  =  2,  3,  .  .  .  ,  k.  Moreover  strict  inequality  holds  for  at  least  one  value 
of  i.  Let  N.  be  the  first  of  these  k  -  1  nodes  to  attain  its  present  label.  Then 
if  D(Nj_j)  denotes  the  label  associated  with  at  the  time  ISL  attained  its 
present  label,  D(N._^)  ¥  D(K_^)  and  thus  D(N^_^)  >  D(N.^)  since  labels  are 
nonincreasing.  Thus 

D(N.)  >  D(Nj_1)  +  e(Nj_1VNJ;F) 

o 

6 

Adding  the  k  1  inequalities,  at  least  one  of  which  is  strict,  and  noting  that 


Ni  -  Nk- 


/■•Tr 

Yd(N.)  >jD(Ni)+^le(Ni_1,Ni;FJ 


O 


But  then  the  cycle  N-,  N^,  .,  ...  .  has  negative  effective  length  relative  to  F. 

O  ;  . 

This  contradicts  lemma  8.  So  no  loop  can  exist,  and  the  proof  is  complete. 


<9  11 


Lemma  11: 


If  D(N)  <  *>for  a  node  N  other  than  A,  an  A/N 
chain  C  exists  whose  links  are  all  approach  links. 


Proof:  Since  N  has  finite  label,  an  approach  link  N2,'N  is  associated 


with  it. :  But  then  has  a  finite  label'knd  an  approach  link  N^,  N^.  One  can 
continue  to  trace  backward  in  this  manner  until  at  some  stage  an  approach  link 


N, 


k’  ^k  i  *s  rea°hed  such  that  -  A.  For  according  to  lemma  10  the  same 


\  node  is  never  reached  twice,  and  there  are  a  finite  number  of  nodes.  Thus 


C  A’  Nk-1’  Nlc-2 


9  •  •  •  9 


N.„,  N;  is  the  desired  chain. 
3  > 


'  Lemma  12: 


/  J 


(( 


Node  A  retains  its  original  label  of  zero.  ^ 


Proof:  Suppose  node  A  is  assigned  a  label  D(A)  =  n  <  0.  0  Let  N,  A  be 
the  approach  link  associated  with  A.  Then  D  (N)  <  °° ,  and  an  A/N.chain  of 


approach  links  A,  ,  N^_  2, 


Ng,  N  can  be  produced  according  to 


lemma  11,  But  then  N,.  A,  N,  , ,  .  .  .  ,  ,N0,  N  is  a  loop,  which  contradicts 

A  k-1  2 

lemma  10.  This  completes  the  proof. 


Lemma  13: 


If  D(N)  <  «*>  the  A/N  chain  C  constructed  in 
Lemma  11  has  the  property  D(N)  >  e(C;F). 
It  follows  that  the  label  associated  with  a 
node  is  never  less  than  the  minimum  ef¬ 
fective  distance  from  A  to  that  node. 
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Proof:  Let  C  be  denoted  by  N^,  N2,  .  .  .  ,  N  where  =  A  and 

))  -  t?  '  n  . 

N  =  N.  Then  for  i  =  2,  3,  , , ,  r  the  following  inequality  holds  according 


to  lemma  9: 


D<N.)  >  D(N._1)  +  e(N._r  N.  ;F) 
Adding  these  inequalities  the  following  is  obtained:. 


^D(N.)  >  Jd(N.)  +  e(C;F) 


°  i=2 


Noting  that  D(N_)  -  0  according  to  lemma  12,  the  desired  result  follows  . 


Now  it  can  be  shown  that  all  labels  remain  constant  after  a  certain 
number  of  repetitions  of  step  b  (see  the  subsection  "Statement  of  the  Algorithm"). 


Lemma  14:  All  node  labels  remain  constant  after  at  most  n  -  1 

iterations  of  (b)  where  n  is  the  number  of  nodes 
in  the  network. 

Proof:  Let  N  be  any  node  other  than  A  such  that  min  [  e(C;F)J  <  «°,  where 

O 

o  - 

the  minimum  is  taken  over  all  A/N  chains  C.  Let  =  N^,  N'2,.  .  .  .  ,  N^.  be 
any  A/N  chain  that  attains  this  minimum.  Clearly  k-<  u  since  no  chain  contains 


the  same  node  twice.  So  C,T  has  at  most  n  -  1  links.  <? 

N 

After  the  first  application  of  step  b  to  all  links,  D(N2)  <  D(N^)  +  e(N^,  N2;F) 

is  obtained  since  N2  is  relabeled  on  the  basis  of  D(N^)  +  e(N^,  N2;F)  unless  an 

til 

even  smaller  label  can  be  assigned.  In  general  on  completion  of  the  i  1  application 


Ji 


of  step  b 


D(N.+i)  <  D(Ni)  +  e(Ni,  Ni+1;F) 


for  i  =  1,  2,  .  .  .  ,  k  -  1,  where  D(N.+1)  is  the  label  of  N.+1  at  the  end  of  the 

til  St  ^ 

i  (or  beginning  of  the  i+1  )  application  of  step  b  and  D(N.)  is  the  label  of  N. 

th 

at  the  beginning  of  the  l  application.  u  Adding  these  k  -  1  inequalities 


k 


k-1 


Td(N.)  <  ^D(N.)  +  e(CN;F) 
1=2  i=l 


Since  D(N^)  =  0  it  follows  that 


•;  / 


D(Nk)  <  e(Cn;F) 


Oil  the  other  hand  D(Nk)  cannot  be  less  than  the  minimum  effective  distance  from 

A  to  N,  since  according  to  lemmas  11  and  13  an  A/N  chain  C  of  approach  links 

0 

can  be  produced  so  that  D(Nk)  >  e(C;F) .  It  follows  that  C  is  an  A/N  cliain  of 
minimum  length,  and  that 


c  D(Nk)  =e(C;F) 

Thus  after  at  most  k-1  applications  of  step  by  D(N)  attains  its  minimum  value. 

C 

l  For  nodes  N  such  that  min  e(C;F)  =°°  ,  D(N)  is  minimized  from  the  start  of  the 
algorithm,  so  that  the  statement  applies  to  all  nodes  of  the  network.  ]  This 


completes  the  proof. 
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The  following  theorem  summarizes  the  foregoing  results;. 


Theorem  4: 


If  the  effective-length  algorithm  is  applied  to 

f  "  ' 

^  .  \>  '.N 

a  network  with  an  ideal  A/B  flow  pattern  F, 

s 

the  process  terminates  (i.  e. ,  all  labels  re- 

'.X 

main  constant)  after  at  most  n  -  1  repetitions 
of  step  b.  Upon  termination  the  final  value 

of  D(N)  for  any  node  N  ~  A  is  equal  to 

((  '  f 

min  [  e(C;F)j,  where  the  minimum  is  taken 
over  all  A/N  chains  C.  Moreover  if  D(N)  <  00 
a  chain  attaining  rnis  minimum  can  be  found 
by  tracing  back  along  approach  links,  as  in 
the  proof  of  lemma  11. 


Thus  the  effective-length  algorithm  possesses  the  properties  that  were 

(, 

j 

assumed  in  the  previous  section.  As  a  by-prOduct  it  determines  the  effective 

distance  from  A  to  every  other  node,  although  only  the  effective  distance  from 

// 

)) 

A  to  B  can  be  used  in  the  procedure.  o 

Note  that  the  sequencing  of  the  algorithm,  whereby  the  links  are  scanned 

in  some  fixed  order  in  step  b,  is  not  claimed  to  be  as  efficient  as  possible. 

6 

Pollack  and  Wiebenson  summarize  a  number  of  efficient  labeling  algorithms 

=.-\  ■  '  (f 

whose  objective  is  to  minimize  1(C)  over  all  A/B  chains  C,  and  it  appears  that 
any  of  these  might  be  adapted  to  minimize  e(C;F).  These  are  generally 
sequenced  in  such  a  way  that  all  links  terminating  in  a  given  node  are  scanned 


together,  and  the  nodes  are  taken  in  a  certain  order.  The  authors  of  this  paper 

.  f/  - 

are  not  in  a  position  to  evaluate  the  relative  efficiencies  of  the  various  schemes 
if  programed  for  a  digital  computer. 
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^  °  ■■  A  , 

GLOSSARY  OF  PRINCIPAL  SYMBOLS  *  * 

‘  ,  f'  . .  c 

’  '  0  " 

/i 

:  .<» 

Symbol 

- 

- 

o  Meaning  '  (  J) 

V/ 

(M,  N) 

Link  joining  nodes  M  and  N  (considered 

& 

\*j 

undirected) 

M,  N 

Directed  link  from  M  to  N 

M^»  Mg,  •  •  •  , 

“k 

An  Ml/Mk  path  if  M^  ~  a  loop  if 

A 

M-  =  M, 

1  k 

c(M,  N) 

Cv  ^ 

Capacity  of  M,  N  (a  nonnegative  integer) 

1(M,  N) 

Length  of  M,  N  (a  nonnegative  real  number) 

(A  1 9  ■^‘2 9  •  •  • 

,Mk, 

o 

Length  of  path  or  loop 

f(M,  N)„  f(N,  M) 

Flow  in  (M,  N)  expressed  relative  to  M,  N 

o 

and  N,  M  respectively 

Y(M) 

Net  output  at  node  M  (flow  output  minus 

input)  < 

V(F) 

Value  of  A/B  flow  pattern  FL  equal  to 

// 

_ 

Y(A)  and  -Y(B)] 

F(+.  G 

\ 

Sum  of  A/B  flow  patterns  F  and  G  v 

f,-..  g 

Difference  between  A/B  flow  patterns 

o 

F  and  G 

C. 

•  •  • 

,Mk>  ( 

The  flow  pattern  obtained  by  assigning  r 

units  of  flow  to  each  Mj,  M.+1  and  zero  , 


to  other  links,,  - 


flow,  summed  oyer  all  links)  c 

o  _ 

Effective  length  of  M,  N  relative  to  flow 

pattern  F  . 


v  Effective  length  of  chain  or  cycle  C 

„ 

relative  to  flow  pattern  F 


// 

length  algorithm. 
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